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THERMAL CAPACITANCE (SLUG) CALORIMETER THEORY INCLUDING 
HEAT LOSSES AND OTHER DECAYING PROCESSES 


T. Mark Hightower,* Ricardo A. Olivares, ^ and Daniel Philippidis* 


ABSTRACT 


A mathematical model, termed the Slug Loss Model, has been developed for describing thennal 
capacitance (slug) calorimeter behavior when heat losses and other decaying processes are not 
negligible, where the temperature time slope decreases noticeably with time. This model is derived 
from heat-transfer and energy-balance principles along with some simplifying assumptions, and 
results in the temperature time slope taking the mathematical form of exponential decay. When data 
are found to fit well to this model, a heat-flux value can be calculated that corrects for the losses and 
may be a better estimate of the cold-wall fully catalytic heat flux, as is desired in arc -jet testing. 

The model was applied to the data from a copper (oxygen-free high-conductivity (OFHC)) slug 
calorimeter inserted during a particularly severe high-heating-rate arc jet run to illustrate its use. The 
Slug Loss Model gave a cold-wall heat flux 15% higher than the value of 2,250 W/cm obtained 
from the conventional approach to processing the data (where no correction is made for losses). For 
comparison, a Finite Element Analysis (FEA) model was created using the commercial software 
program COMSOL Multiphysics, where conduction heat losses from the slug were simulated. Very 
close agreement between the temperature -versus-time result of the FEA model and the actual data 
was obtained by a unique solution of both the heat flux and heat-loss coefficient specified in the 
FEA model. A sensitivity analysis showed that this value of heat flux as determined from the FEA 
model was detenninable to within ±1%, and was also found to be in close agreement with the heat 
flux detennined by the Slug Loss Model. The FEA model accounted for temperature-dependent 
physical properties, heat capacity with the Shomate equation, and thermal conductivity with a simple 
linear fit. Although this work was applied to arc -jet test applications of slug calorimeters, it should 
have general applicability wherever slug calorimeters are used and losses are significant. The 
idealized theory of slug calorimeters for no losses and constant physical properties is also developed 
as a necessary basis for the general derivation presented. 


Engineer, Thermophysics Facilities Branch, M/S 229-4, Ames Research Center, Moffett Field, CA. 

1 Engineer, Thermophysics Facilities Branch, M/S 213-8, Ames Research Center, Moffett Field, CA. 

* Masters Student, Department of Mechanical and Aerospace Engineering, San Jose State University, San Jose, CA. 
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I. INTRODUCTION 


A slug calorimeter determines heat flux by measuring the rate at which a slug of material heats up 
while subjected to a heat source. Arc jets are ground test facilities that are used to produce heating 
and flow environments similar to those experienced in planetary atmospheric entry in order to test 
spacecraft thermal protection materials and systems. Slug calorimeters are used for calibration of 
arc -jet test conditions (ref. 1). For arc -jet applications the slug is usually made of oxygen-free high- 
conductivity (OFHC) copper. Figure 1 shows an assembly drawing of an arc -jet slug calorimeter. 

Heat losses from the slug to its holder are of concern, and are difficult to control under high-heat- 
flux conditions (ref. 2). The American Society for Testing and Materials (ASTM) standard (ref. 1) 
recommends using the slug cool-down slope after removal from the heat source as an indication of 
the losses during the heating phase. This approach assumes that the conditions causing the losses are 
the same whether subjected to or removed from the heat source — an assumption that evidence 
suggests is not always valid. The ASTM standard also recommends that these losses be less than 5% 
of the heating rate in order to be neglected. 

This current work explores how losses during the heating phase can be determined more directly 
based on the behavior of the temperature time slope during the heating phase, where the temperature 
time slope decreases more rapidly than can be explained by an increase in the heat capacity of the 
copper slug alone. Furthermore, this work explores how these losses can be accounted for in the 
computed value of heat flux. Correcting slug calorimeter results for heat losses to achieve more 
accuracy has been recommended in the literature (ref. 3). 


Body 



Set Screw 


slug i 




Ruby Sphere 


Thermocouple 

Wire 


Holder 


Figure 1. Assembly drawing of hemispherical slug calorimeter. 
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To provide the necessary groundwork, Section II covers slug calorimeter theory where no heat 
losses and constant physical properties are assumed. Then Section III covers slug calorimeter theory 
where losses and variable heat capacity with temperature are taken into account. Section IV shows in 
detail the application of the theory of Section III to real data from one arc -jet run, and Section V 
shows and compares results from a FEA model of the same run. Finally, Section VI presents 
conclusions. 


II. NO HEAT LOSSES AND CONSTANT PHYSICAL PROPERTIES 


A study in which this specific problem was fully derived was not found, although one or more must 
surely exist. Instead a textbook was found wherein the applicable method of solution was clearly 
explained, allowing for application of the method to this problem (ref. 4). 

Consider a slug calorimeter, a right circular cylinder of solid OFHC copper, with diameter D and 
length L, initially at constant uniform temperature T 0 . At t = 0 unifonn constant heat flux q is applied 
at one end, the front face, while the back face and cylindrical surface are idealized to be perfectly 
insulated, as shown in figure 2. 

Therefore, the heat flow can be modeled as one-dimensional unsteady-state heat transfer. The 
temperature of the slug is measured at the back face with a thermocouple. Define x = 0 as the front 
face and x = L as the back face. 


The one-dimensional unsteady-state heat-balance differential equation is 

d 2 T _ pc p dT 
dx 2 k dt 

The thermal diffusivity is defined as 

k 

a = 

P°p 


(2.1) 


( 2 . 2 ) 



Figure 2. Boundary conditions for idealized model. 
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The boundary conditions are 


dr(0,t) _ q 
dx k 

dT(L,t) _ Q 

dx 


The initial condition is 


T(x,0) = T o 


(2.3) 


(2.4) 


The main steps in the solution of this boundary-value problem by the method of separation of 
variables are outlined as follows. All physical properties are assumed to be constant. 


First express the solution as the sum of two parts, the steady-state solution and the transient solution. 
Overall solution = steady-state solution + transient solution 

(2.5) 

T(x,t') = v ss (x,t') + w(x,t) 

The “steady-state solution” here is the solution that applies when t becomes large, where the temper- 
ature everywhere in the slug is increasing at the same rate, such that 

dT 

— = constant (2.6) 

dt 

The transient solution is a correction to the steady-state solution for when t is small, and it 
approaches zero as t becomes large. Through derivation the steady-state solution was determined to 
be 


v ss (x,t) = T a + 


qt qL qx~ qx 

Lpc 3 k 2 Lk k 


(2.7) 


That this solution is correct can be verified by noting that it satisfies: the differential equation (2.1), 
the steady-state criterion of equation (2.6), and equation (2.8) based on an overall energy balance on 
the slug from t = 0 to t = t. In words this equation states that the heat added to the slug must be 
reflected in its average temperature. 


T =To + - 


qt 


Lpc L 


1 ^ 

= — J Tdx 


( 2 . 8 ) 
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Substituting v ss (x,t) for T in the integral of equation (2.8) and evaluating the integral gives 


Y l y 1 

-\v ss (xj)dx = -\ 
Li Li 


f . r 2 'N 

T [ qt { qL | qx qx 


Li 


o 

Lf 




Lpc 3 k 2 Lk k 


dx = 


T+- 


qt 




Lpc 


f , h 

t + ^L 

v l /*pj 
^ qt ' 


T+- 


V 


Z/7C 


l 

+ — 
L 

+ 0 


1 L ( T 2 

1 r qx qx 

l{ , 3k 2 Lk k 


dx ■ 


f (ik+ 

2 

qx 

\ 

qx 


2 Lk 

k V 


dx = 


p J 


(2.9) 


Applying the change of variable from T to w to the differential equation gives 


d 2 w 1 dw 


dx 2 a dt 


( 2 . 10 ) 


The boundary conditions in w are 


dw(0, t) 
dx 

dw(L,t ) 
dx 


= 0 
= 0 


(2.11) 


The initial condition is 


»<*,0 ) = z£l + £_^ 

2Lk k 2k 


(2.12) 


It is the change of variables that allows w to be solved by the separation of variables method, 
because the boundary conditions in w are homogeneous. So the separation-of-variables method 
proposes that 

w(x, t ) = </)(x)d(t) (2.13) 
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Substituting into the partial differential equation produces two ordinary differential equations: 


d w 

-— = < /,”(x)0(t) 
ox 

dw 

—=<p(x)m 

at 

</>"(x)0(t) = -</)(x)0Xt) 
a 

rw _ 1 0\t) _ 12 

<p(x) a 0{t) 

0"+A>=O 
0'+oc; l 2 0 = O 

The negative sign in front of A 2 is needed because without it a trivial solution is 
general solutions to these differential equations are 

(f>(x) = a cos Ax + b sin Ax 
0(t) = e a2}t 

The boundary conditions become 

ox 

h«M =W )<w=o 

dx 

dw(Lj ) =</>'( L )0(t) = Q 
dx 

Since 0(t ) = 0 would give a trivial solution, the boundary conditions become 

^’( 0 ) = 0 
<P'(L) = 0 

Apply the first boundary condition 

^’(0) = 0 

(f) '(x) = -a A sin Ax + bA cos Ax 
^'(0) = -<7/lsin(0) + h/lcos(0) = 0 
b = 0 

/. (j)(x) = a cos Ax 


(2.14) 


obtained. The 


(2.15) 


(2.16) 


(2.17) 


(2.18) 
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Apply the second boundary condition 


(j)\L) = -aAsmAL = 0 
0 nn 

A/n — 

L 

n = 1, 2, 3... 

At this point the solution to the partial differential equation in w can be expressed as 

oo 

w(x, t) = ^ a n cos 


^ UlTtX ^ 


n = 1 


tin 

-a\ — \t 


V L J 


Now it remains to apply the initial condition in order to solve for the constants a n . 

/ nnx 


w(x,0) = ^<3„cos 


n = 1 


V L J 


1 

1 

N> 

qx_ 

qL 

J 2 Lk 

k 

~2k ~ 


Equation (2.21) is a problem in Fourier series, which is solved by 

2 qL 1 


2 

a n =-JgO)cos 


^ Mix '' 


0 


V L j 


dx = -~ 


kn 2 n 2 


The transient solution then becomes 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


(2.22) 


. 2 qL 1 

M<x - ,)= ^S^ cos 


nnx 

V L j 


(2.23) 


Finally, combining equations (2.7) and (2.23) yields the overall solution to the partial differential 
equation: 


T(x,t) = 


T | qt [ qL | qx_ qx 
Lpc 3 k 2 Lk k 


2 qL -A 1 

t/ — rcos 

kn 2 n 


f nnx'' 


\ L j 


Factoring out common terms and applying equation (2.2) gives 


T(x,t) = T 0 + ^ 

k 


at 1 1 f x ^ 


- H 1 — 

V 3 2 




«=1 


r _2 Lu 2 

L TC n= \ Tl 


Z~ cos 


( nnx\ _Q (y^ * 


\ L j 


(2.24) 


(2.25) 
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As a check on equation (2.25), it is worth comparing it to the equivalent solution given in the classic 
text by Carslaw and Jaeger (ref. 5). What makes comparing the solutions problematic is that their 
text defines the x coordinate as zero at the back face of the slug and L at the front face. So call the 
coordinate in the Carslaw and Jaeger text x c to distinguish it from the x of this paper, and make a 
change of variable in equation (2.25) according to the following equation: 

x = L-x c (2.26) 


Making this change of variable and simplifying yields the following equation, which is exactly 
equivalent to that given in the Carslaw and Jaeger text when T 0 is set equal to zero (ref. 5). 


TM-T.+ f §44 f 



(2.31) 



The average temperature of the slug T ave is given by 


T = T (?) — T + 

1 ave 1 ave ' L ' 1 o 


qat 
kL L 


1 1 

= — J Tdx 


(2.32) 


The following difference equations can be written: 


T _ T = }_qk^_4qL 


I- 


1 


Jl2n^X t 


2 k kff (2n - 1) 

j _ 1 qL | 2qL yi (-1)" ^ p L 
b 6 k k7T 2 n 


(2.33) 


The summations in all of these equations approach zero as t becomes large. 

It can be shown that if enough time has elapsed for the heat to have just penetrated to the backside of 
the slug, the summation terms are well represented by truncating the summations with the first term, 
for t greater than or equal to the penetration time t p . 


qat \qL 2qL 
b 0 kL 6 k kn 2 

Take the derivative and set it equal to zero to solve for t p : 


(2.34) 


dT b _ qa 2 qa - I •- 

dt kL kL 


(2.35) 


*p = 


L 2 


ajr 


-ln(2) 


Now do a more general treatment for response time ?«. Make a distinction between q, npu t to the front 
face and q indicated based on the slug calorimeter equation applied at the back face. 


q = q, 


input 


dT h 

Q indicated P^p ^ 


(2.36) 


Vindicated 

( T dT h > 

Lpc p 2 

\ dt j 

_ pc pOC 

2 pc p a 

V input 

V input 

k 

k 


=\-2e^ ^ 
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f 


Ir 


L 2 


-In 


cck~ 


Q indicated 


1- 


9, 


( 2 . 37 ) 


input J 


For practical purposes, the response time calculated when ^ md,caled =0.99 should be sufficient 


q, 


input 


elapsed time for the heat-flux determination from the back-face temperature to begin to be valid. 


L 2 


f 


In 


1-0.99 


( 2 . 38 ) 


In the next section this value tRo .99 is used as the response time necessary to ensure that the slug has 
reached steady state, where the temperature profile has established its parabolic shape. 


III. HEAT LOSSES AND VARIABLE HEAT CAPACITY - SLUG LOSS MODEL 


In this section the theory is developed to handle the case of slug calorimeters where heat losses from 
slug to surroundings (as well as other possible decaying processes) are taken into account. This 
theory is developed in order to create a model to fit to real data, where losses are often evident by a 
T-vcrsus-/ slope that decreases with time more than can be accounted for by the increasing heat 
capacity of the copper slug with temperature alone. 


A heat balance on the slug with losses gives 


input - output {i.e. losses ) = accumulation 


qA _(J^lA = Mc “T 


R 


dt 


( 3 . 1 ) 


The surroundings (slug holder) to which the slug loses heat is assumed to be at the same temperature 
as the initial temperature of the slug, T a , and this surrounding temperature is assumed to remain 
constant. R[ a is defined as the apparent resistance to heat loss, which is coupled with c po , the constant 
value of c p at temperature T 0 . These parameters are coupled so that equation (3.1) can be analytically 
integrated. Later the actual resistance to heat loss, R/, will be solved for when heat capacity as a 
function of temperature, c p (T ave ), is introduced. At this point the effect of increasing c p with 
temperature is accounted for as an additional apparent loss incorporated into the value of R/ a . 
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It is desired to get equation (3.1) in terms of T b , the temperature that is recorded as data. Assume that 
enough time has elapsed so that the parabolic temperature profile through the slug has been 
established, a condition that should be the case for t > tRo.99 from equation (2.38). From equa- 
tion (2.33) for t > tRo.99, the summation drops out as zero and gives 

T„-l \=^r ( 3 . 2 ) 

6k 

Solving this equation for T ave gives 

t„=i ;+ f ^ < 3 . 3 ) 

6k 


Substituting this result into equation (3.1) gives 


A 

q 

l {Me, 


L 

6kR Ia Mc po j 


+ - 


R ,a MC poJ 


KMc po 


dT ; 

dt 


( 3 . 4 ) 


Define two constants: 


a = 


A 


L 


V y MC po 6 kR la MC poJ 


+ - 


R la Mc po ) 


R, Me 

la po 


( 3 . 5 ) 


Now equation (3.4) can be written as: 


a ~bT h 


dT l 

dt 


( 3 . 6 ) 


This equation integrates to 


T - — 

Hi fit , 


e ~ b (t-h) + bL 


( 3 . 7 ) 


This equation can be used to fit actual slug calorimeter data recorded during an arc -jet run by taking 
the following steps: 

First determine the time when the slug arrived at the position where the heat flux is to be measured. 
Then add tRo.99 from equation (2.38) to this time to get ti, the initial time for which the data will be 
processed. This step ensures that steady state has been reached by time t\. Next, detennine fc, the 
final time for which the data will be processed, by taking the time just before the slug begins to be 
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removed from the heat source. Now there is a range of 7) r vcrsus-/ data points from tj to t 2 to fit with 
equation (3.7). (If the slug has been left exposed to the heat source longer than is necessary, an 
alternate approach might also be used to select t 2 . In such a case, adding tm .99 or perhaps 1.5 times 
t R( i 99 to ti to get t 2 might be a better approach. This approach would help ensure the validity of the 
assumption of constant surrounding temperature T 0 up to the value of t 2 where the model is being 
applied.) 

Next assume a value of “b” and plot 7), versus e H ' h 1 , which should be linear according to the 
equation. Next determine the value of “b” that gives the best linear lit. This determination can be 
made by solving for “b” that maximizes the R" value (coefficient of determination) for the lit. Now 
that “b” is known, “a” can be solved for from the intercept of the lit. Finally, Tbif,t can be solved for 
from the slope of the lit. The model equation fit to the data is now complete and ready to be 
compared to the data. 

After the data are processed and “o” and “b” have been detennined, q can be solved for by solving 
equation (3.5) for q: 


a — bT 


Mc po ( a ~bT 0 ) 


f 



A - 

L 

1 


yMCpo 6kR la Mc po j 


A 


1 — 


L 


V 6kR la A J 


An alternate but equivalent equation for q is 

Me. 


q 


A 


a -h 

r T _ 

V 

l ° 6k ), 


( 3 . 8 ) 


( 3 . 9 ) 


Take equation (3.7) and go back in time to t 0 , the theoretical time when the step function of “q on” 
would have occurred had it been a perfect step function. This process goes back in time while 
keeping the temperature profile of the slug at steady-state. It is an idealization that differs from the 
actual behavior, but whether idealized or actual, the energy that heats up the slug from t 0 to ti is the 
same. At t = t 0 , T ave = T 0 . According to equation (3.2) 

TM = T 0 ~ ( 3 - 10 ) 


When going back in time while maintaining a steady-state parabolic temperature profile, the 
temperature at the back face at t 0 has to be less than T 0 so that T ave ~ T 0 . Substituting equation (3.10) 
into equation (3.7) gives 


T M = T.~ = 

6k 


1 SI fit 


A 


V 


e ~b(t 0 -h) + Q_ 

b 


( 3 . 11 ) 
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Solving this equation for t Q gives 


1, 

t o= t i-T ln 

b 


f T qL 

0 6k b 


a 


1 bl fit 


(3.12) 


Equation (3.7) gives Tb as an analytical function of time, which can be written as Tb(t) to make the t 
functionality explicit. By using the linear relationship between T ave and Tb of equation (3.3), T ave can 
now also be expressed as an analytical function of time, T ave (t). 


(f 


T ave (t) = 


vv 


T 


e - b (t-h) + ^_ 


+ 


qL 

6k 


The derivative of T ave (t) with respect to time can also be taken analytically: 


dT ave {t) 

dt 


■—b 


l blfit 


i) 


(3.13) 


(3.14) 


Now c p can be introduced as a function of T ave into the differential equation (3.1) as well as the 
actual loss resistance, Ri. 


qA _ T a ) = Mc ^ (T ^ (0) dT ave (. t ) 


R, 


dt 


(3.15) 


This equation can now be solved for R/ as a function of time: 

( T m (t)-T a ) 


R,(t) : 


f dT (t) 

qA-Mc(T ave (t)) a 


dt 


J 


(3.16) 


It remains to express c p for copper as a function of T ave by the Shomate equation. 

c p (T„,) = A + BT „ + CTJ + DTJ + JL 

± ave 

where 

A = 2.789933xl0 2 — B = 4.42 1789x1 0" 1 J 


kgK 


kgK 


C = -4.918 152x1 CT 4 ^— D = 2.19879x1 0" 7 J 


E = 1.079706x10 


kgK 
eJK_ 
kg 


kgK 4 


(3.17) 
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The Shomate equation is available at the National Institute of Standards and Technology (NIST) 
web site (ref. 6). The equation given by NIST was modified to get the temperature units as shown in 
equation (3.17), and to get it on a mass basis rather than a mole basis, the conversion of the atomic 
weight for copper, 63.546 g/gmole, was applied. 

So now all the elements that allow the calculation of R/(t) by equation (3.16) have been assembled. 

Other useful equations can also be written. Two equations can be written for heat flux based on the 
temperature-time slope and heat capacity, one based on Tb and the other based on T ave . 


qslopeTh ’ A dt 

(3.18) 

l slope 1 ave v ' 

(3.19) 


Equation (3.18) is essentially equivalent to the current ASTM method where no correction for losses 
is attempted if time is chosen to be early in the steady-state region, i.e., t = tj. Equation (3.19) is 
similar, but is probably slightly more accurate, being based on T ave rather than Tb. Keep in mind that 
in the above equations where temperatures and derivatives of temperatures are explicitly expressed 
as functions of time, these functions of time all stem from the fit of the model to the data, 
equation (3.7). 

By comparing equation (3.19) to equation (3.9), noting equations (3.10) and (3.6), recognizing that 
the time derivatives of T ave and 77, are equal by inspection of equation (3.3), and that Tave(to ) - To , it 
can be shown that 


9 d slopefavc \^o f 


Me., 


A 


a—b 


T - 


1L 'll WlMj 


6k 


JJ 


A 


dt 


(3.20) 


Mc r J dT l At.)\_Mc rc dT a ,,(t„) 


A 


dt 


J 


A 


dt 


Finally, the following equations can be written for heat- flux loss and fraction loss as functions of 
time: 


d,aJ t ) = d-d s lopeTaJ t ) = d 
FracLoss(t) = 1 - p 


M c p (T ave (t )) dTjn 
A dt 

Mc p (TJt)) d T ave (t) 


qA 


dt 


(3.21) 


(3.22) 
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Although the method developed in this section is based on the assumption of heat losses, resulting in 
exponential decay of the temperature time slope, the possibility that other phenomena might mimic 
exponential decay should not be dismissed. For example, the front surface of the slug might start out 
being fully catalytic, but become less so very rapidly upon insertion into the hot oxidizing plasma 
flow as to result in a decaying net heat flux to the slug. Under such a scenario the Slug Loss Model 
developed here would still provide some correction for these “losses.” The next section presents the 
application of this model for including losses to the data from one run. 


IV. ANALYSIS OF DATA FROM RUN IHF187R025 


Figure 3 shows a plot of the data, where the slug back-face temperature and stagnation pressure are 
plotted as functions of time. The stagnation pressure is used in this run to determine when the slug 
reached centerline and when it left centerline, in order to be able to select the data to be processed 
according to the method. 

Careful analysis of the stagnation pressure data shown in the figure indicates that t c in and t c out are 
325.992 s and 327.102 s, respectively, where t c in is the time when the slug has just reached 
centerline and t c out is the time just before the slug starts to leave the centerline. T 0 is determined from 
these data to be 302.35 K. 



Figure 3. Back-face temperature and stagnation pressure versus time. 
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The physical properties of the copper slug are: 

Density at 298 K, assumed a constant throughout calculations: 

p = 8,925.7^- (4.1) 

m 

Thermal conductivity at 298 K, assumed a constant throughout calculations: 

W 

k = 385.2 (4.2) 

mK 

Heat capacity calculated by the Shomate equation (3.17) at T a : 


c =385.615 

" kgK 


(4.3) 


Other data for the slug follow: 


D = 0.0078 Im (4.4) 

M = 0.004529 kg (4.5) 

A = 0.25 nD 1 = 0.000047906 m 2 (4.6) 

L= — = 0.010592 m (4.7) 

pA 


Combining equations (2.2) and (2.38), setting c p = c po , and substituting for these defined values 
gives 


R0.99 


pc po L 2 

kn 1 


In 


f 2 ) 

v 1-0.99 j 


= 0.538 s 


(4.8) 


Add the value of tRo .99 to t c », to get tj = 326.53 s, which defines the lower end of the time span for 
which the data will be processed. The closest value to this in the data turns out to be ti = 326.532 s. 
The upper end of the time span, C, is simply equal to t c out = 327.102 s. With ti and C defined, the 
data to process are shown in table 1 . The data-acquisition sampling rate was about 67 data points per 
second. 
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TABLE 1. BACK-FACE TEMPERATURE VERSUS TIME DATA FROM ti TO t 2 


t, S 

Tb, K 

t, s 

Tb, K 

t, s 

Tb, K 

t, s 

Tb, K 

326.532 

660.7955 

326.682 

744.1884 

326.832 

825.295 

326.983 

903.1487 

326.547 

668.8717 

326.697 

752.7801 

326.847 

833.1806 

326.997 

909.7945 

326.562 

677.0599 

326.712 

760.8678 

326.862 

841.5302 

327.012 

917.442 

326.577 

686.4049 

326.727 

768.6463 

326.877 

848.8967 

327.027 

925.1961 

326.592 

694.5348 

326.742 

777.3664 

326.892 

856.1385 

327.042 

932.1872 

326.608 

703.5074 

326.757 

785.3169 

326.907 

864.8431 

327.057 

939.6873 

326 622 

711.217 

326.773 

794.0334 

326.922 

872.3267 

327.072 

947.399 

326.637 

719.5663 

326.787 

801 .3908 

326.937 

879.3349 

327.087 

954.5518 

326.652 

728.1516 

326.802 

809.1195 

326.952 

887.2804 

327.102 

961 .6053 

326.667 

736.1388 

326.818 

818.2387 

326.967 

895.0791 




Fitting this data to equation (3.7) gives the best linear fit to T b versus e h(t when b = 0.29160 s' 1 , 
where the R“ value of the fit is maximized at 0.99999. The Solver function in an Excel spreadsheet 
was used to solve for b. 


When b is determined, “a” can be calculated according to the equation 

a = intercept (r h vs ) b = 766.76 


K_ 

s 


( 4 . 9 ) 


And T bl f lt can be calculated according to the equation: 


T u f , =slope(r, Vi ) + fU 660.32 K 


( 4 . 10 ) 


Ri a can be calculated according to the equation: 


R„ 


bMc, 


1.964 — 
W 


( 4 . 11 ) 


And q can be calculated according to equation (3.8): 


Me (a-bT) W W 

q = — ^ ^ = 26, 005, 000 — = 2, 600 — I 


A 


1 - 


L 


6kR, a Aj 


in 


cm 


( 4 . 12 ) 


This value is about 15% higher than the value of 2,250 W/cm reported by the facility test engineers, 
where losses were not taken into account. 
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And t 0 can be calculated according to equation (3.12): 


/ =t, — In 


( T qL 

0 6k b 


a 


bl fit 


= 325.792 5 


( 4 . 13 ) 


The fictitious Tb(t 0 ) of equation (3.10) can be calculated: 

TM = T„-^- = mK ( 4 . 14 ) 

6k 

The plot in figure 4 shows how the model fits the data. 

Notice how the fit is extremely good, with error of fit to data of ± 0.07%. Also note that the fit is 
clearly not linear; this fact is easily seen by holding a straight edge to it. Figure 5 shows a plot of the 
losses versus time as determined from equations (3.21) and (3.22). 

Notice how at ti, the earliest point at which the slope would be eligible to be used with the current 
ASTM method (without corrections for losses), the losses are greater than 10%. And the losses only 
get larger with time. 

For comparison, in the test engineer’s report for this run, a heat loss of 345.4 W/cm“ was reported 
during the cooling phase, after the slug was removed from the flow. 

Figure 6 shows a plot of the actual loss resistance, R/, versus t, calculated according to equa- 
tion (3.16). 



Figure 4. Back-face temperature - fit compared to data. 
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— q losses 

— FracLoss 


Figure 5. Losses (per cm 2 slug frontal area) versus time. 



t, s 


Figure 6. Actual loss resistance versus time. 
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V. FINITE ELEMENT ANALYSIS MODEL 


A Finite Element Analysis (FEA) model was created using the commercial software program 
COMSOL Multiphysics, COMSOL, Inc, Burlington, Massachusetts, USA. The slug was modeled 
using three-dimensional (3-D) tetrahedral elements. Figure 7 shows the mesh. 

The heat-loss routes from slug to holder for this model were defined as follows (reference to figure 1 
is helpful for understanding this description): The actual slug is isolated from its holder by six ruby 
spheres of 1.5875-mm diameter, three arranged equally spaced around a concentric circle at the back 
side of the slug, and three arranged equally spaced around the circumference of the slug about 
midway from front to back. Based on measurements of pronounced dimpling observed on a slug, 
from the ruby spheres being forced into the copper surface of the slug during use, it was estimated 
that the area of contact between each ruby sphere and the slug could be represented as a circular area 
with diameter of 0.6 mm. 

Therefore, for the FEA model, the surface areas of the slug that make contact with the ruby spheres 
were modeled as 0.6-mm-diameter surface regions, with a heat-transfer coefficient of h rejecting 
heat to a temperature of T a , the initial temperature of the slug and assumed constant temperature of 
its holder. These features of the model are shown in figure 8. Note that for the purposes of this 
model possible heat losses through the air gap between slug and holder via convection and radiation 
were not considered. 



Figure 7. Model mesh. 
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Figure 8. Boundary conditions showing the locations of the entering and exiting heat fluxes 

in the model. 


The slug has the material properties of copper. The heat capacity varies with temperature according 
to the Shomate equation (3.17). Thermal conductivity (ref. 7) was fit to a linear equation (5.1). 

k = c x T + c 2 

q =—0.071098— ^ (5.1) 

mK 

W 

c 2 =422.915 

mK 

A smoothed Heaviside function (flc2hs from the COMSOL Library) was used to create ramp-up and 
ramp-down times for the heat flux applied to the front face of the model. The total simulation time 
of the model is 3 seconds, with a 0.01 -second time step. Figure 9 shows a temperature fringe plot at 
the end of a run (at time t = 3 seconds). 

Runs were made with the model for various values of q input, h, and duration of heat pulse, until a 
very close agreement with the actual data was obtained with a value of q input = 2,600 W/cm 2 . Fig- 
ure 10 shows Tb versus t for both this COMSOL solution and the actual data, where the curves lie 
almost exactly on top of one another, and figure 1 1 compares the q values as calculated from the 
slug calorimeter equation applied to the back face, i.e., per equation (3.18) as applied numerically to 
the data. 
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Max: 1045.76 K 



I 1040 K 
1030 K 
1020 K 
1010 K 
- 1000 K 
- 990 K 
1980K 


- 970 K 


1 940 K 
Min: 937.511 K 


Figure 9. Temperature fringe plot of the model at time t = 3 seconds. 



Figure 10. 7^ vs. f-COMSOL model and actual data compared. 
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The best fit of the model to the data is obtained with a unique combination of q and h. A sensitivity 
analysis showed that the values of q and h so detennined are kn own to about ±1% and ± 10%, 
respectively. The temperatures of the front and back faces were plotted as a function of time; they 
are shown in figure 12. 



Figure 11 . q versus t - COMSOL model and actual data compared. 


Temperature vs Time 



— Front Face 
— Back Face 


Figure 12. Temperature vs. time plots for the center of the front and back faces. 
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Alternate COMSOL models were used that were axisymmetric (2-D) where the ruby sphere contacts 
were modeled as thin concentric rings of contact area. These models gave essentially the same 
results as the 3-D model already presented, but had the advantage of running much faster. In the data 
presented in the following discussion, the COMSOL models are all based on the 2-D models, but the 
results would be expected to be nearly identical if the more time-consuming equivalent 3-D models 
were run. 

COMSOL model runs were done for the following four cases: no losses and constant physical 
properties, no losses and variable physical properties, losses and constant physical properties, and 
losses and variable physical properties (the 2-D version of the 3-D case already presented). In 
addition, the ideal partial differential equation solution of equation (2.30) was calculated. These five 
cases were all based on a q input value of 2,600 W/cm . Finally, the Slug Loss Model fit of equa- 
tion (3.7) was calculated based upon the COMSOL data for the “losses and variable physical 
properties” case. Figure 13 compares all six of these cases. 



Ideal PDE Soln Perf Step 

No Loss Const Phys Props 

No Loss Var Phys Props 

— Loss Const Phys Props 

Loss Var Phys Props 

Slug Loss Model 


Figure 13. Results of all models compared. 
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The ideal partial differential equation solution has a perfect step function, so its initial transient 
solution is quicker than the “no losses and constant physical properties” COMSOL case that has a 
smoothed step function, but as expected, these two solutions parallel each other beyond the transient. 

The “no losses and variable physical properties” case shows more deviation from the ideal case than 
the “losses and constant physical properties” case, showing that variable heat capacity contributes 
more to the dropping slope during the heating phase than the losses do. 

The Slug Loss Model solution fit the “losses and variable physical properties” COMSOL case over 
the steady-state region, from t = 0.76 s to 1.3 s, the region over which the lengthwise temperature 
profde is roughly parabolic. As expected, the Slug Loss Model gives a value of Tt(t 0 ) less than T 0 
according to equation (3.10), where t 0 calculates to be about 0.04 s according to equation (3.12). 

When the Slug Loss Model was run on the data from the “losses and variable physical properties” 
COMSOL case, it gave a q of 2,465 W/cm 2 , about 5% less than the 2,600 W/cm 2 value that was 
determined from both the Slug Loss Model applied to the actual data and the COMSOL lit to the 
actual data. The reason for this difference appears to be evident in figure 1 1 . The downward slope of 
the actual data is steeper over the time interval where the Slug Loss Model was fit than the 
downward slope of the COMSOL model fit over the same interval. The steeper downward slope is 
indicative of more losses, which is inferred as a higher q by the Slug Loss Model. Thus, the Slug 
Loss Model appears to have a bias towards underpredicting q (i.e., not fully correcting for losses) if 
the data were to have perfectly followed the assumptions of the COMSOL model. Based on this 
finding, it is best not to apply the Slug Loss Model blindly without also looking at the data plotted in 
the form of q versus t as in figure 11 to be sure that the data show the basic characteristics upon 
which the Slug Loss Model is based. 

Note that in the actual slug calorimeter data case under analysis in this paper, the front face of the 
slug actually melted at the end of the run. This result did not appear to have a discernable impact on 
the data. No attempt was made to incorporate this melting in any of the models. 


VI. CONCLUSIONS 


A mathematical model for determining the heat flux with a slug calorimeter with high heat losses, 
tenned the Slug Loss Model, has been presented. For a particular high-heat-flux arc -jet run, the Slug 
Loss Model detennined a heat flux that is 15% higher than what has been calculated using the 
conventional approach, ft was found to be in good agreement with the FEA model, reinforcing the 
idea that losses are occurring. This Slug Loss Model is intended to be an additional tool in the 
thennal analyst’s toolbox that can be applied for high losses. 
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